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We derive constraints on Wilson coefficients of dimension-six effective operators 
probing the b — > s transition, using recent improved measurements of the rare 
decays B s — > [i + fj,~, B — > K/j> + fi~ and B — > K*fi + fi~ and including all relevant 
observables in inclusive and exclusive decays. We consider operators present in 
the SM as well as their chirality-flipped counterparts and scalar operators. We find 
good agreement with the SM expectations. Compared to the situation before winter 
2012, we find significantly more stringent constraints on the chirality-flipped coeffi- 
cients due to complementary constraints from B — > K[i + and B — > K*fi + fi~ and 
due to the LHCb measurement of the angular observable S3 in the latter decay. We 
also list the full set of observables sensitive to new physics in the low recoil region 
of B K*n + fi-. 

1. Introduction 

Rare B meson decays probing the flavour-changing neutral current b — > s transition are sen- 
sitive to physics beyond the Standard Model (SM). Independently of a concrete model, the 
effect of new degrees of freedom that are sufficiently heavy compared to the decaying B mesons 
can be described by the modification of Wilson coefficients of local operators in an effective 
Hamiltonian of the form 



Since the set of relevant operators and the dependence on the Wilson coefficients can be dif- 
ferent for the various experimental observables probing this Hamiltonian, a combined analysis 
of all available experimental constraints is mandatory to obtain meaningful bounds on the in- 
dividual coefficients and to determine the room left for new physics (NP). Since the startup of 
the LHC, such constraints have become significantly more stringent, thanks to improved mea- 
surements of exclusive decays. In the last few months, the LHCb and BaBar collaboration have 
presented updated analyses of branching ratios and angular observables in B — > [i + , 
while the LHCb, ATLAS and CMS experiment have presented strong bounds on the branch- 
ing ratio of the very rare B s — > decay. The purpose of this paper is to update the 
model-independent constraints on Wilson coefficients first presented in IT]. As in the previous 
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analysis, we go beyond comparable recent studies [2j|7| by considering the SM operator basis 
as well as their chirality- flipped counterparts and by including the most complete set of observ- 
ables in inclusive and exclusive decays. With respect to [l], we emphasize the most important 
improvements: 

• We take into account the new data on B — > fi + fj,~ and B s — > released by 
LHCb and BaBar in 2012. 

• We include the decays B — > Kfi + fi~ and B s — > n + fi~ in the determination of Wilson 
coefficients constraints. 

• We list the full set of B — > K* /i + angular observables sensitive to NP at high dilepton 
invariant mass q 2 , pointing out that there are 5 observables only sensitive to right-handed 
currents. This is particularly relevant for S3 recently measured by LHCb and CDF, which 
we now include in our analysis also at high q 2 . 

• We discuss the impact of the direct CP asymmetry in B — > X s ^f on the constraints, 
considering the sizable theory uncertainties due to long-distance contributions. 

Next, we will discuss the relevant observables and their dependence on NP in sec. [2j focusing 
on the improvements relative to [I] . The resulting constraints on individual Wilson coefficients 
are presented in sec. [3j Sec. [4] is devoted to a discussion of the potential size of NP effects in 
observables yet to be measured, which is useful to assess the prospects of future measurements. 
We summarize our results in sec. [5] Finally, appendix |A| contains a detailed discussion of our 
method to average measurements from different experiments and compares these averages to 
our theoretical predictions. 



2. Observables 

We consider NP contributions to the operators 

O (0 = ^{sa^Pn^F^, of = 9 -^{sa, u T a P R(L) b)G^\ 

Of = {s^P L {R)b){hn) , Og = (s7,P L{R) b)(h^) , 

Of = ^(sP R{L)b )(U) , 0<? = —(SP R{L) b)(£ l5 £) , (2) 

and define all the Wilson coefficients at a matching scale of 160 GeV. In the SM, the primed 
coefficients as well as Cs,p are negligibly small. For numerical values of the Wilson coefficients 
we refer to the appendix of (lj. Our normalization choice for the scalar and pseudoscalar 
operators is chosen to make their Wilson coefficients dimensionless and scale invariant. We 
note that, by factoring out the CKM elements in the definition of ([!]), a non-zero imaginary 
part for a Wilson coefficient corresponds to CP violation beyond the CKM phase. 

We consider all relevant observables in inclusive and exclusive B meson decays that are 
sensitive to the Wilson coefficients of the operators in Q. In particular, in addition to the 
observables considered in [I] - the branching ratio of the inclusive B — > X s "f decay, the time 
dependent CP asymmetry in the exclusive B — > if* 7 decay, the branching ratio of the inclusive 
B —7- X s £ + l~ decay, the branching ratio and angular observables in the exclusive B — > K* /i + /i~ 
decay and the branching ratio of the rare B s — > {i + [i~ decay - we also include the branching 
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ratio of the exclusive B — > Kfj, + fi~ decay, that provides valuable complementary information on 
the right-handed Wilson coefficients, and the direct CP asymmetry in B — > Xgj, Acp(b — > s"f). 

Our treatment of BR(B — >• X s ^f) and B — > K*^f exactly follows our analysis in [Tj. We 
slightly refine our treatment of B — > X s l + l~ by adding a 10% theory uncertainty due to 
the cut on the invariant mass of the hadronic final state [§]. In view of the current large 
experimental uncertainty, the impact on our conclusions is marginal. Improvements in our 
treatment of B — > K* /U + ^~ and B s — > fi + n~ on both the theory and experimental side, as 
well as our analyses of B — > Kp^ and Acp{b — > sj), are described in detail below. 

2.1. B s 

Recently, it has been pointed ou t [28[[29] that the large width difference in the B s system, 
y s = rs s Ar s /2 = (8.8 ± 1.4)% [ 30 1 , has to be taken into account when comparing theory 
predictions of BR(i3 s — > with the experimental results. The corresponding correction 

is not universal but depends on possible NP contributions to B s — > In presence of NP 

one has 

BR( B ,^+,-) ^Vl 2 + |P| 2 )xfl + ,, R ^>- R p f 2 ' V(l + fe ). (3) 



BR(B S -> /x+/i-) S M V ) V \S\ 2 + \P 

where 



q _ m Bs (C S - C's) 
2m, CfM 



p = m Ba (C P - C'p) | (Cio - C[ ) 



m 



2 

B a I ^10 °10 



(4) 



In (|3j), BR(B S -> /i+/i~)sM and BR(B S -» refer to the SM and NP predictions for 

the branching ratio extracted from an untagged rate that can directly be compared to exper- 



iment 1 28 29 . Using an average of the most recent lattice determinations of the B s meson 
decay constant ||l}|33|, f Bs = (227 ± 4)MeV §4], we fincQ 



BR(B S M /Osm = (3.32 ± 0.17) x 10~ y . (5) 



The experimental sensitivities start to close in on the SM value. While CDF [16] finds a 
small excess in B s — > /x + [i~ candidates and gives a two-sided limit on the BR(S S — > fi + fi~) 
at 95% C.L., LHCb |13|, CMS [14] and ATLAS [15] quote at 95% C.L. only upper bounds 
on the branching ratio. At the 2cr level all results are consistent and we perform a naive 
combination using the available A% 2 distribution for BR(1? S — > /x + /i~) from CDF as well as 
the CL S distributions from ATLAS, CMS and LHCb. We obtain 

BR(B S = (2.4 ± 1.6) x 10- 9 , (6) 

which is in good agreement with the SM expectation. While this clearly leads to strong con- 
straints on models with potentially large contributions to the scalar or pseudoscalar operators, 
contributions to the Wilson coefficients C\q, C[ were known even before the measurement ^ 
to be unable to enhance the branching ratio above about 5.6 x 10 -9 due to constraints from 
other rare semi-leptonic decays, as shown in [l]. As a consequence, the experimental precision 
for this decay is only now starting to become sensitive to models modifying the branching ratio 
by means of the semi-leptonic operators Oio and 0' w . 

1 Neglecting the y s correction we find BR(£? S — > (i + h~)sm = (3.05 ± 0.15) x 10~ 9 , which corresponds to the 
CP averaged branching ratio. 
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Observable Experiment 


SM prediction 


10 4 xBR(£ -»• X s j) 3.55 ±0.26 < 
A CP (b^s~f) (-1.2 ±2.8)% ! 

Q,^ fl 1 fi + fl 99 ( 
iJK*^ — U.J.O it U.ZZ ; 

10 9 x BR{B S -»• /x + At~) 2.4 ±1.6 


} 
) 
j 

L3- 


16 


3.15 ±0.23 
(1.0 ±1.6)% 
(-2.3 ± 1.6)% 
3.32 ±0.17 


10 
11 
12 




K , , N [1,61 1.63 ±0.50 

10 6 xBR(5 -> x^+r 

0/1^-1-019 


17 


18 


1.59 ±0.19 
0.23 ±0.07 


8, 
20 


19 


[1,6] 1.28 ±0.12 
10 7 xBR(£ -»• ff*+r) [14.18,16] 0.41 ±0.06 

[1 R 99 Ql n /IQ -1- D f)7 

[io, zz.yj u.^yztu.ui 


21 


- 24 


1.29 ±0.30 
0.43 ±0.10 
0.86 ±0.20 


[1,6] 1.90 ±0.19 
10 7 xBR(£ -> K*l+£-) [14.18,16] 1.12 ±0.11 

M fi 1 Ql 1 97 -I- n 1 9 


21 


- 25 


2.49 ± 0.64 
1.13 ±0.41 
1.34 ±0.62 


[1,6] 0.63 ±0.06 
(F L )(B ^ K*£+£-) [14.18,16] 0.34 ± 0.06 

[16 191 32 ± 06 


99 


25 


0.77 ±0.04 
-27 0.37 ±0.20 
0.34 ±0.25 


[1,6] 0.09 ±0.05 
(^FB>(^ -> K*£ + £-) [14.18,16] -0.46 ±0.05 

[16,19] -0.39 ±0.06 


22 


25 


0.03 ± 0.02 
-27 -0.41 ±0.13 
-0.35 ±0.13 


[1,6] 0.05 ±0.08 
(S 3 )(B ^ K*£+£-) [14.18,16] 0.02 ± 0.08 

[16,19] -0.21 ±0.08 


25 


26 


(-0.3 ± 1.1) 10~ 2 
-0.14 ±0.08 
-0.22 ±0.10 


[1,6] 0.09 ±0.39 
(A 9 )(B -> K*£ + £-) [14.18,16] 0.18 ±0.25 

[16,19] -0.20 ±0.34 


26 


(1.5 ±2.4) 10" 4 
(-0.7 ± 1.7) 10~ 5 
(-0.5±1.2)10" 5 



Table 1: Experimental averages and SM predictions for the observables used in the fit. SM 
predictions that lack a reference are based on our calculations. 
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Obs. 


low q 2 


high q 2 


BR 


CV 9 in 1 C7 Q 10 


Cq in, Cq in 


F L 


C*7 Q . C7 n in 

1 j" ' ( ,y, iu 


r" 

w 1 n 

y,iu 


s 3 


CJ 

in 
1 ,1U 


r" 

y , ±u 




CV 10 , Ci 1 n 


r" 

wn in 

y , ±u 


s 5 


C7 Q , O7 in 


Cq , Cq 1 n 

17 3 y , iu 


A F b 




C9,10; Cg,io 


A 7 


CV,10> ^7,10 




As 


CV,9, C79 10 




A 9 


u 7,10 


u 9,10 



Table 2: Observables in the angular distribution of B — > K* and Wilson coefficients they 

are most sensitive to in the low and high q 2 regions. 



2.2. B K*n+n 



The B — > K*fj, + fi decay offers a multitude of observables that are sensitive to NP effects 
(see e.g. [2,35-37]). Updated measurements of the branching ratio and angular observables 



in B — > K*[i + fi have recently been presented by BaBar 23 , 27 and LHCb 



24 



25 



In 

addition to the branching ratio, forward-backward asymmetry j4fb and the K* longitudinal 
polarization fraction Fl, LHCb also strongly improved the constraints on the CP-averaged 



angular observable S3 first constrained by CDF 1 26 
to NP in right-handed currents. Compared to [l 



, which is tiny in the SM but highly sensitive 
, we now also include the £3 constraint at 



high q and find it to give a relevant constraint on the primed Wilson coefficients. We also 
include the angular CP asymmetry Ag, that is very sensitive to CP violation in right-handed 
currents, both at low and high q 2 . Currently, only weak bounds on this observable exist from 
CDF 26 and their impact on our results is small. The impact of a future measurement of Ag 



by LHCb will be discussed at the end of section |4j 

Concerning the high q 2 region, we also emphasize that with our operator basis, eq. ([2]), 
there are more independent observables at high q 2 than in the case of the SM operator basis. 
For example, while it has been shown in j2] that Fl is "short-distance free" at high q 2 in 
the case of the SM operator basis and can be used to extract form factor ratios from the 
data |38| , this is no longer true in the presence of NP contributions to C' 7 9 10 . 



Consequently, 

Fl at high q 2 can be considered a probe of right-handed currents. The same is true for the 
observable £4 defined in |36| and for the CP asymmetry ^8- The CP asymmetry A7 is instead 
"short-distance free" at high q 2 even in the presence of chirality-flipped coefficients, while S5 
at high q 2 is sensitive to NP even with SM-operators only. The sensitivity to NP in the various 
observables is summarized for convenience in table [2j 

On the experimental side, both LHCb and BaBar observe some hints for a non-zero isospin 
asymmetry. However, since the theoretical expectation for the isospin asymmetry is rather 
limited compared to the experimental uncertainties even in the presence of NP (39 j, we feel 
justified in neglecting it in the following and use both the B + and B° data in our averages. 
Details on the averaging procedure and the comparison of experiment vs. theory are discussed 
in appendix \A\ 
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Finally, we note two minor numerical refinements of the SM predictions compared to [l] : 

• We use an updated value for the neutral B meson decay constant fg = 190(4) MeV [34] 
relevant for the normalization of B — > K* form factors, the effect being an increased SM 
prediction for the branching ratio in the low q 2 bin (but small with respect to the theory 
uncertainty); 

• For the extrapolation of form factors to high q 2 , where we use the results of [40] , we now 
take into account the error correlations between the fitted expansion parameters. This 
leads to a very small increase in the estimated theory uncertainties. 

Our predictions for all the relevant observables at low and high q 2 are listed in table [l] As 
we take into account error correlations for the B — > K* form factors in the low q 2 region |36| , 
some of our theory predictions are considerably more precise than the values in w. 



2.3. B Kn+n~ 

The branching ratio of the B —> K/j, + fi~~ decay has been measured by the BaBar, Belle, CDF 



and LHCb collaborations [21-24 
has recently been analyzed in 



Its impact on constraints on the Wilson coefficients CV^io 
6j. Here, we include the BR(i? — > Kfi + fi~) in our analysis, 
discussing also the chirality-flipped Wilson coefficients 9 10 . Since the B — >■ K transition 
does not receive contributions from an axial vector current, the primed Wilson coefficients 
enter the B — > Kfi + /j,~ observables always in conjunction with their unprimed counterparts as 
(Ci + C-). This is complementary to the B — > K*fj, + /j,~ decay and is useful to constrain the 
chirality-flipped operators. Neglecting lepton mass effects and (pseudo)scalar operators, the 
differential branching ratio can be written as [41] 



dBR r B G 2 F a 2 e \V th V t %\ 2 A 3 / 2 2 2 

W = 2%5 m 3 3- (1^1 +\F V \), (7) 



where 



F A = (C 10 + C' 10 )f + (q 2 ) , Fy = (Cf + C' 9 )f + (q 2 ) + —^—{ C f + C' 7 )f T (q 2 ) , (8) 

A = m% + m% + q 4 - 2(m 2 B m 2 K + m 2 B q 2 + m 2 K q 2 ) . (9) 

In ^ the Wilson coefficients are evaluated at the scale of the b quark, fi ~ mj. The effective 
Wilson coefficients are defined for instance in |36|. In addition, we take into account 
non-factorizable 0(a s ) corrections proportional to form factors [42]. For the B — )• K form 
factors, we use the results of (40] , which are based on a fit to a LCSR calculation valid at 
low q 2 and lattice results valid at high q 2 . We take into account all the error correlations 
among the fit parameters given in an unpublished updated version of |40] 2 [ In view of the 
sizable resulting form factor uncertainties, we neglect the impact of possible higher-order non- 
perturbative corrections, which are estimated to be at the few percent level [43[|44] . Similarly 
to the case of B — > K*fi + fi~, we ignore isospin breaking effects, which are expected to be well 
below the experimental sensitivity, and consequently use the average of B° and B + data for 
the experimental branching ratios. Our resulting SM predictions for the branching ratio in the 
different q 2 bins are shown in table [l] 

2 We thank Aoife Bharucha for providing us with the correct covariance matrices. 
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2.4. A CP (b sj) 



The direct CP asymmetry in the B — > X s j decay probes CP phases in the Wilson coefficients 
C-j and C' 7 of the magnetic operators. It arises first at NLO and therefore depends also in a 
non-trivial way on the Wilson coefficients Cg and C§ of the chromo-magnetic operators. This 
is in contrast to all the other observables sensitive to ci that we consider in this work that 
all depend on the same combination of ci and Cg . 

The experimental world average for Acp(b — > S7) is dominated by data from Belle [45] and 



BaBar 46 and reads [9] 

A C p{b ^ sj) cxp = (-1.2 ±2.8)% , (10) 

which is consistent with CP conservation. 

The SM prediction for Acp(b — > 57) suffers from large hadronic uncertainties [TT] 



- 0.6% < A C p(b -> s 7 ) exp < 2.8% , (11) 

which makes it difficult to identify small NP contributions in this observable. Nonetheless, 
Acp(b —> 57) can be used to constrain large CP violating NP contributions which would 
otherwise be allowed by measurements of CP-averaged quantities like branching ratios. 

In our numerical analysis, we use the expressions given in [II] that we extend to include 
also the right-handed Wilson coefficients C' 7 and Cg. To estimate the theoretical uncertainty 
in presence of NP, we consider the uncertainties in the so-called resolved photon contributions 
treating the model estimates for the A parameters as Gaussian la ranges, as well 



given m 
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as the uncertainties in the remaining contributions coming from scale variation and m c /mb 
given in [47]. We add all uncertainties in quadrature. Using this procedure, we obtain the SM 
prediction quoted in table [T] 



3. Model-independent constraints on Wilson coefficients 

Combining all the experimental averages and theory predictions of the observables listed in 
table [TJ we construct a x function and use it to obtain model-independent constraints on the 
Wilson coefficients. Our statistical treatment follows [I], with one minor refinement: while 
previously treating all theory uncertainties as uncorrelated, we now treat the uncertainties of 
the two adjacent high-g 2 bins for all observables in B — > K* / u + and B — > Kfi + fj,~ as 100% 
correlated, which we found to be approximately fulfilled for generic values of Wilson coefficients 
and which we checked to lead to slightly more conservative (looser) constraints than assuming 
them to be uncorrelated. 

We first discuss the constraints on single complex coefficients or on pairs of real Wilson 
coefficients. Figures [l] and [2] shows contours of A\ 2 = 4 (~ 2a constraints) from individual 
processes and A\ 2 = 1,4 (~ 1,2a constraints) from the combined constraints. We make the 
following observations. 

• At the 95% C.L., all Wilson coefficients are compatible with their SM values. For the 
total x 2 f or the SM values of the Wilson coefficients we find xV^dof = 21.8/24. This 
value improves only slightly in presence of NP. 

• For the coefficients present in the SM, i.e. C7, Cg and C10, as well as for C Q and C' 10 , the 
constraints on the imaginary part are looser than on the real part. 
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Figure 1: Individual 2a constraints on the unprimed Wilson coefficients from B — > X s £ + l~ 
(brown), BR(B ->■ X s j) (yellow), A C p{b — > sj) (orange), B — > K*j (purple), B ->■ 
K*/j, + fi~ (green), B — > Kfi + n~ (blue) and B s — > fi + fi~ (gray) as well as combined 1 
and 2a constraints (red). 
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Figure 3: Constraints on the scalar Wilson coefficients cjp and C^ p at the 1 and 2a level from 
BR(i? s — > fi + fi~), assuming no new physics in C^g. 

• For the Wilson coefficients cjg, the new constraint on B s — > fi + pT is starting to become 
competitive with the constraints from B — > K^*' p + pT . 

• The constraints on C 9 and C[ from B — > Kp + fi~ and B — > K* /i + fi~ are complementary 
and lead to a more constrained region, and better agreement with the SM, than with 
B —7- K* /i + pT alone. 

• A second allowed region in the C7-C7 plane characterized by large positive contributions 
to both coefficients, which was found to be allowed e.g. in [lpj, is now disfavoured at 95% 
C.L. by the new B — > K* p + fj,~ data, in particular the forward-backward asymmetry. 

The Acp(b — > S7) constraint limits the size of allowed imaginary contributions to Cj. Its 
impact on C' 7 is marginal, so we refrain from showing it. 



The second point above can be understood from the fact that in the branching ratios and CP 
averaged angular observables giving the strongest constraints, only NP contributions aligned 
in phase with the SM can interfere with the SM contributions. As a consequence, NP with 
non-standard CP violation is in fact constrained more weakly than NP where CP violation 
stems only from the CKM phase. This highlights the need for improved measurements of CP 
asymmetries directly sensitive to non-standard phases, such as the T-odd CP asymmetries in 
B -> K*n+n~. 

The new experimental bounds on BR(i? s — > p + can be used to constrain the combinations 
of scalar and pseudoscalar Wilson coefficients Cs,p — C' s p . The resulting bounds are shown 

in figure 3, assuming c[q to be SM-like. The right plot showing the constraint for Cs — 
C' s = Cp — C'p can be easily reinterpreted in the frequently studied framework of the minimal 
flavor violating MSSM in the large tan/3 regime where Cs — —Cp and C' s ~ C' P ~ C10 — 
C' w ~ 0. Complementary bounds on the combinations Cs,p + C' s p can be obtained from the 
B —7- K[i r [T' decay, as has been done recently [7j. However, given the strong bound from 
B s — > /U + /i _ , this is only relevant if one considers Cs,p + C' s P 3> Cs,p — C' s P . 

The constraints on the individual Wilson coefficients can be translated into bounds on the 
suppressing scale of the flavour- violating dimension-six operators. Table [3] shows the con- 
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Operator A [TeV] for \a\ = 1 |c,| for A = 1 TeV 



+ — +i —i + - +i —i 



O t = 


?f(sa^P R b)F^ 


69 


270 


43 


38 


1.6 


10" 


-4 


9.7 


10" 


-6 


4.2 


10" 


-4 


5.3 


10" 


-4 


0' 7 = 




4G 


70 


78 


47 


3.6 


10" 


-4 


1.6 


10" 


-4 


1.3 


10" 


-4 


3.5 


10" 


-4 


o 9 = 




29 


G4 


21 


22 


1.2 


10" 


-3 


2.4 


10" 


-4 


2.2 


10" 


-3 


2.0 


10" 


-3 


o' 9 = 


{si^p R b){hn) 


51 


22 


21 


23 


3.8 


10" 


-4 


2.1 


10" 


-3 


2.2 


10" 


-3 


1.9 


10" 


-3 


o w = 


(s^P L b)(£^ 5 £) 


43 


33 


23 


23 


5.4 


10" 


-4 


9.2 


10" 


-4 


1.9 


10" 


-3 


1.9 


10" 


-3 


O' 10 = 


(s^p R b)(ir-f 5 e) 


25 


89 


24 


23 


1.7 


10" 


-3 


1.3 


10" 


-4 


1.7 


10" 


-3 


1.9 


io- 


-3 


— 


^t s {sP R[L) b){U) 


93 


93 


98 


98 


1.1 


10" 


-4 


1.1 


10" 


-4 


1.1 


10" 


-4 


1.1 


10" 


-4 


o P = 


^(sP R b)(I l5 £) 


173 


58 


93 


93 


3.3 


10" 


-5 


3.0 


10" 


-4 


1.1 


10" 


-4 


1.1 


io- 


-4 


o' P = 




58 


173 


93 


93 


3.0 


10" 


-4 


3.3 


10" 


-5 


1.1 


10" 


-4 


1.1 


io- 


-4 



Table 3: Lower bounds (at 95% C.L.) on the NP scale related to the relevant dimension six 
operators, assuming the coefficient c% in the effective Hamiltonian T-L e g = — aOi/A 2 
to be +1, —1, +i or —i as well as upper bounds on the coefficients Cj for a NP scale 
of A = 1 TeV, depending on Cj/|cj|. 



straints on the scales Aj in the effective Hamiltonian^] 



j=7,9,10,5,P ^ ' 

assuming |cj| = 1, as well as the constraints on the q, assuming A = 1 TeV. The constraints 
were obtained by varying only one coefficient at a time and correspond to A% 2 = 4. In both 
cases, we show the constraints for Cj/|cj| = +1, or —i. 

On average, the obtained constraints are slightly weaker compared to the constraints that one 



finds for the suppressing scales of dimension-six AF = 2 operators that lead to B s mixing 48 
Exeptions are the magnetic operator if it interferes with the SM destructively in B — > X s "f or 
scalar operators if they interefere constructively in B s — > fi^ pT . Those operators are probed 
to scales above 100 TeV, as are operators that induce B s mixing. 



4. Allowed effects and future prospects 

The constraints in the previous section were derived assuming only one Wilson coefficient at a 
time, or the real parts of two coefficients, to deviate simultaneously from the SM. In the generic 
case, where all Wilson coefficients are allowed to deviate from the SM, cancellations may occur 
which render some of the constraints ineffective. On the other hand, even if one takes into 
account such cancellations, the current data already put indirect limits on observables which 
have not been measured to a good precision yet. 

To obtain such predictions for allowed regions in the presence of generic NP, we perform 
a Bayesian Markov Chain Monte Carlo (MCMC) analysis of six different scenarios where 

3 The constraints on the d can be easily translated to constraints on the d defined in eq. [l]by dividing them 
by (lTeV) 2 i|f lW sI f^2 « 8.2 x 10" 4 . 
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Re(C 10 ) 



Figure 4: Constraint on C' 10 at 68% (green dashed) and 95 % C.L. (green solid) from a hypo- 
thetical measurement of Ag with the SM central value and the experimental errors 
equal to the current LHCb errors on Sg, compared to the current global constraint 
on C 10 (red). 

subsets of the Wilson coefficients or all of them have been varied, assuming flat priors for the 
real and imaginary parts of the Wilson coefficients, and using e~ x I 2 as likelihood, with the x 2 
function described above^](for more details see [l]). We then determine the Bayesian posterior 
credibility regions for each observable (for central values of the theory parameters) in each 
scenario. These ranges can be interpreted as allowed regions in the different scenarios (with 
or without right-handed currents, with or without CPV beyond the CKM phase) and give an 
indication of the prospects of future measurements of these observables. 

In table [4], these fit predictions are shown at 95% C.L. for the branching ratio of B s — > 
and for various angular observables in B — > K*fi + fj>~ sensitive to NP. We observe in particular 
that 

• a suppression of BR(i? s — > fj, + fj,~) below 10~ 9 requires NP in both left- and right-handed 
currents, i.e. in both C\$ and C{ J^] 

• The T-odd CP asymmetries ^7,8,9 can still exceed 10% (for Aj even 30%) in the presence 
of non-standard CPV. For effects in Ag at high q 2 and Ag, right-handed currents are 
required. 

While A§ and Ag at high q 2 are null tests of the SM and therefore promising inspite of 
sizable relative uncertainties, we point out that S3 at high q 2 is non-zero even in the SM and 
afflicted with a large uncertainty (cf. table [I]), so that the large variation of the central value 
in the presence of right-handed currents shown in the last line of table [4] is spoiled by the badly 
known SM value. This would change if the relevant form factors could be estimated more 
precisely, e.g. by means of lattice QCD. 

Finally, we want to emphasize the importance of measuring the angular CP asymmetry Ag 
in B — > K*fi + fi~. As is well known, Ag is highly sensitive to NP in right-handed currents 

4 In obtaining the predictions in table Ul we have not included Aqp (b — > S7) in the \ 2 function, since C ( 7 ' ] and 
Cg enter this observable in a different linear combination compared to all other observables, so in a generic 
fit of Wilson coefficients, the constraint could always be compensated by adjusting Cj\. 

5 We remind the reader that here we assume (pseudo-) scalar currents to be absent. 
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< 0.28 





< 0.22 





< 0.35 


\(Aa)[i,o]\ 





< 0.14 





< 0.15 





< 0.21 


K^s) [14.18,16] 1 











< 0.13 





< 0.12 


K^9>[l,6]l 











< 0.09 





< 0.13 


[14.18,16] I 











< 0.22 





< 0.20 


(£3) [1,6] 








[-0.02,0.08] 


[-0.04,0.09] 


[-0.06,0.15] 


[-0.07,0.13] 


(S3) [14.18,16] 


-0.14 


-0.14 


[-0.16,0.01] 


[-0.17,0.03] 


[-0.18,0.01] 


[-0.17,0] 



Table 4: Predictions at 95% C.L. for the branching ratio of B s — > and predictions for 

angular observables in B — > K*[i + n~ (neglecting tiny SM effects below the percent 
level). The columns correspond to 6 scenarios with real (1R) or complex (C) new 
physics contributions to the operators CV^io and/or C' 7 91 q. 



at low q 2 , and we stressed above that this is also true for the high q 2 region (cf. table [2]). 
Moreover, Ag can be obtained from a merely one-dimensional angular distribution. Using the 
conventions of [36], the distribution in the angle eft reads 



d(r + f) Mr + f) 1 



d(j) dq 2 



dq 2 



— [1 + S 3 cos(2<j 

Z7T 



' + Ag sin(2</>)] , 



(13) 



that is, Ag can be extracted from an untagged sample 35 just as S3. Indeed, in 2011 CDF 
presented the first measurement of Ag (26) , denoted Ai m by CDF. The preliminary results for 
the quantity also denoted Aj m recently presented by LHCb 25 correspond insteacj^Jto the CP- 
averaged angular observable Sg, whose sensitivity to NP is very limited 36 . To demonstrate 
the constraining power of Ag on CP-violating right-handed currents, we show in figure [4] the 
constraints from a hypothetical measurement of Ag with SM central value (i.e. zero), and the 
experimental errors equal to the current errors of LHCb's measurement of Sg. As the figure 
shows, such measurement, which should be possible even with the current LHCb dataset, would 
significantly reduce the allowed size of NP in Im(CJ ). 



5. Conclusions 

Using the recent improved measurements of rare decays probing the b — > s transition, we have 
derived updated constraints on Wilson coefficients of dimension-six effective operators. We 
went beyond comparable recent studies by including the most complete set of experimental 
observables and NP operators, including also the chirality-flipped counterparts of the operators 
present in the SM. Compared to the predecessor study [l], in addition to including the new 

6 Thomas Blake and Nicola Serra, private communication. 
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experimental data, we made several improvements, the most significant of which being the 
inclusion of the decays B — > Kfi + fj,~ and B s — > fJ, + fjT as well the direct CP asymmetry in 
b —7- 57 in our constraints. 

The most important results can be summarized as follows. 

• At the 95% C.L., the Wilson coefficients are still compatible with their SM values. 

• The branching ratio of B — > Kfi + ' \i~ constitutes a particularly important constraint on 
the chirality- flipped Wilson coefficients C' 9 10 , for which it is complementary to B — > 

• The new LHCb measurement of the angular observable S3 in B — > K*n + fi~ constitutes 
a significant constraint on right-handed currents. It is sensitive to NP also in the high 
q 2 region. 

• A measurement of the CP asymmetry Ag is expected to provide a significant constraint 
on Im(Cg 10 ) even with the current LHCb dataset, if measured to be compatible with the 
SM. 

• In the high-g 2 region of B — > K* n + [x~ , there are five angular observables only sensitive 
to right-handed currents: Fl, S3, S4, As and Ag. 

In the near future, improved constraints - or chances to uncover physics beyond the SM! 
- will be facilitated by improved measurements of the branching ratios of B — > Kfi + /j,~, 
B — > X s £ + £~ and B s — > and the remaining CP-symmetric and asymmetric angular 

observables in B — > K*fi + fj,~. As we have also seen, several observables in exclusive decays 
are already dominated by theory uncertainties, so progress in lattice calculations of hadronic 
form factors would be very welcome. 
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A. Averaging procedure and theory vs. experiment in 
B K<* W 

In this appendix, we describe our method to obtain the averages of experimental measure- 
ments of B — > K* ^ and B — > K[i + fi~ observables quoted in table [TJ Given a number of 
measurements of the same experimental quantity (like a partial branching fraction), we first 
symmetrize asymmetric statistical and/or systematic errors using the prescription of ref. [49]. 
Then, we perform a weighted average of the symmetrized individual results. We use a slightly 
modified procedure in the case where experiments quote separate results for the charged and 



neutral B — > K^fi + fj,~ decays. As mentioned in section 2.2, we assume the isospin asymmetry 
to be small with respect to the experimental uncertainties in all observables. However, in some 
cases, in particular the BaBar B — > K*fi + fi~ data at low q 2 , the results for the charged and 
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Figure 5: Symmetrized experimental results and averages obtained using the procedure de- 
scribed in appendix [A] for the B — V K*fi + fi~ branching ratio, Fl and Afb in the 
three q 2 regions considered. The green bands are our theory predictions with la un- 
certainties. In the cases where we perform the average of B° and B + data ourselves 
(for details on the procedure see text), we also show the separate results as dashed 
error bars. 
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Figure 6: Symmetrized experimental results and averages obtained using the procedure de- 
scribed in appendix |A] for the B — > K[i + n~ branching ratio in the three q 2 regions 
considered. The green bands are our theory predictions with la uncertainties. 



neutral modes show a significant difference. In these cases, we first perform the average of the 
charged and neutral mode results of a single experiment using the PDG averaging method, 
i.e. rescaling the error by a factor of y/x 2 , and use this rescaled error in the weighted average 
with the other experiments. If the charged and neutral mode results agree within experimen- 
tal uncertainties or if the experimental sensitivity to one mode is much better than to the 
other, we use the combination provided in the experimental publications, if available. Our SM 
predictions for the B -4 K^fi + fi^ branching ratios in table [l] have been obtained using the 
neutral B lifetime, so in the combination of experimental results, we have rescaled the charged 
B branching ratios, if necessary. 

Figure [5] confronts the individual (appropriately symmetrized) experimental measurements 
of the branching ratio, Fl and ^4fb in B — > K*fj, + fj,~ and our average to the SM prediction. 
We highlight the following points, 

• The theory uncertainties are particularly large for the branching ratio and the high-g 2 
angular observables. In these cases, the uncertainty is already dominated by theory 
and progress in constraining NP will only be possible with better control over the form 
factors, which might come from lattice QCD. 

• In the low-q 2 angular observables, there is some tension between BaBar and LHCb re- 
sults, in addition to the tension among charged and neutral B results observed by both 
experiments. 

Figure [6] shows the same comparison for the B — > K/j, + fj,~ branching ratio. Also here, the 
uncertainties are dominated by theory and require progress on the lattice calculation of form 
factors. We also observe a tension in the highest q 2 bin at a level of somewhat less than 2a, 
which however does not have a big impact on the constraints since the adjacent bin agrees well 
with the SM expectation and NP typically affects both bins simultaneously. 
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